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Abstract 

With Monte Carlo simulations, we systematically investigate the depinning phase transition in 
the two-dimensional driven random-field clock model. Based on the short-time dynamic approach, 
we determine the transition field and critical exponents. The results show that the critical expo- 
nents vary with the form of the random-field distribution and the strength of the random fields, and 
the roughening dynamics of the domain interface belongs to the new subclass with Q ^ Q oc ^ ( s 
and Qioc ^ 1. More importantly, we find that the transition field and critical exponents change 
with the initial orientations of the magnetization of the two ordered domains. 
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I. INTRODUCTION 



In recent years the domain-wall dynamics in disordered media has attracted much at- 
tention, which is important from both experimental and theoretical perspectives. A crucial 
feature of the domain-wall motion in disordered media is the depinning phase transition 
at zero temperature. A typical theoretical approach to the domain-wall dynamics is the 
Edwards- Wilkinson equation with quenched disorder (QEW) [l-6]. It is usually believed 
that different variants of the model belong to a same universality class. 

The QEW equation is a phenomenological model, without microscopic structures and 
interactions of the materials. In two dimensions, the domain wall is assumed to be a single- 
valued elastic string. An alternative, possibly more realistic approach, is to construct lattice 
models at the microscopic level. The simplest example is the two-dimensional (2D) driven 
random- field Ising model (DRFIM) 7Hll| . In the DRFIM model, overhangs and islands are 
created in the dynamic evolution, and the domain wall is not single-valued. Very recently it 
has been demonstrated that the DRFIM model does not belong to the universality class of 



the QEW equation, due to the dynamic effect of overhangs and islands 
are closer to experiments 
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- Il2j , and its results 



16j ]. The DRFIM model does not suffer from the theoretical 



self- inconsistence as in the QEW equation, and may go beyond the QEW equation, for 



example, to describe the relaxation-to-creep transition 
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13| . To distinguish the DRFIM 



model and QEW equation, one needs accurate measurements of the critical exponents. In 



this respect, the short-time dynamic approach has been proven to be efficient 
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Traditionally domain walls in magnetic materials are viewed as Ising-like [22J. However, 
this assumption has been challenged by the complex Bloch walls or Neel walls. In recent 
years many activities have been devoted to the orientation dynamics of the magnetization 
distribution, in experimental, theoretical and numerical studies 23M271]. Particular attention 
is also drawn to the depinning transition {28 1. Recently domain walls converging at different 
angles are observed in the thin films of BiFeC>3, PbTiC>3 and PbZro.2Tio.sO3 [29|, l30 |. 

Up to now, most theoretical studies of the depinning phase transition of the domain in- 
terface in magnetic materials are concentrated on the Ising-like systems, such as the DRFIM 
model. In both theories and experiments, it is desirable to go beyond. For example, one may 
consider the XY model, in which the spins are planar unit vectors. However, the pinning 
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phenomenon does not exist in the 2D XY model. Possibly, this might be connected to the 
fact that there are no real ordered states at finite temperatures in the 2D XY model. 
The p-state clock model is a discrete version of the XY model 
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32| . In the equilibrium 



state, the 2-state clock model is just the Ising model; for p — 4, the clock model is decoupled 
into two Ising models; for p = 3, it is equivalent to the three-state Potts model; for p — > oo, 
the clock model approaches the XY model 



. For 4 < p < oo, there is a quasiliquid 



intermediate phase between two transition temperatures T\ and T 2 (T\ < T 2 ). Here T\ is 
the transition temperature from the ordered phase to the quasiliquid phase, and T 2 is that 
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from the quasiliquid phase to the disordered phase 

With Monte Carlo simulations, and based on the short-time dynamic approach, the 
purpose of this paper is to study the depinning phase transition in the 2D p-state clock 
model with quenched random fields, emphasizing the dynamic effect of the orientations of 
the magnetization of the two ordered domains. Numerical results are mainly presented for 
p = 6, and also extended from p = 2 to 12. In Sec. II, the model and scaling analysis are 
described. In Sec. Ill, numerical results are presented. Sec. IV includes the conclusions. 



II. MODEL AND SCALING ANALYSIS 



A. Model 



The 2D p-state driven random-field clock model (DRFCM) is defined on a square lattice 
by the Hamiltonian 

u = - J E & ■ Sj - J2(h + ■ ^ 

<ij> i 

= -J cos(9i - 9 j) - H ^2 cos 9i -^2 hi cos(ipi - 9i), (1) 

<ij> i i 

where the spin Si = (cos^j, sin 6^) = (S itX , S itV ) is a planar unit vector at site i of a rectangle 
L x x L y , and 6i is the angle associated with the orientation of the spin, taking discrete 
values 9i = 2irqi/p with g« = 0, . . . , p — 1. The sum < ij > is over the nearest neighbors, 
if is a homogeneous external driving field in the x direction, and tpi is the angle associated 
with the orientation of the random- field hi, taking discrete values (pi = 2imi/p. Here rij 
is a random integer within an interval [0,p — 1]. The random field hi may be distributed 
in different forms. A typical example is a uniform distribution within an interval [—A, A], 
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3- 0, Q, 



following Refs. [7|, |10|, We take A from 1.0J to 2.5J and set J = 1. Another example 



is a Gaussian distribution with standard deviation a and zero mean. The form of the 
random-field distribution and the strength of the random fields may change the order of the 
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depinning phase transitions and also its universality class 

We consider a dynamic process with a semiordered initial state with a perfect domain 
wall in the y direction. As time evolves, the domain wall propagates and roughens, therefore 
it looks like an interface. We first concentrate on the uniform distribution of the random 
fields, and then extend the calculations to the Gaussian distribution. Our simulations are 
performed at zero temperature with lattice sizes L y = 256, 512, 1024, and 2048 up to t max = 
5120. Total samples for average are about 20000 for different random fields. Main results 
are presented with L y = 1024 and simulations of different L y confirm that finite-size effects 
are already negligibly small. L x is taken to be sufficiently large so that the boundary is not 
reached. Typically, L x = 300 to 500. Statistical errors are estimated by dividing the total 
samples in a few subgroups. If the fluctuation in the time direction is comparable with or 
larger than the statistical error, it will be taken into account. 

For an even p, spins in the initial state take the orientation 0l = on the sublattice at the 
left side, while Or = ti on the sublattice at the right side. For an odd p, spins in the initial 
state take the orientation 0l = on the sublattice at the left side, while Or = 7i(p — l)/p 
on the sublattice at the right side. Further, we also change the initial orientation Or on the 
sublattice at the right side, to examine its influence on the depinning phase transition. A 
periodic boundary condition is used in the y direction, while spins at the boundary in the x 
direction are fixed. To eliminate the pinning effect irrelevant for the disorder, we rotate the 
lattice such that the initial domain wall orients in the (11) direction of the lattice [?], lo| . 

After preparing the initial state, we randomly select a spin and a new state % and flip 
it if the total ener gy d ecreases after flipping, and at least one of the nearest neighbors is 



already flipped [10 
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31] . In other words, only spins at the interface are allowed to be flipped 



A Monte Carlo time step is defined as L x x L y single spin flips. In Fig. (TJ snapshots 
of the spin configurations are displayed for the 6-state DRFCM model with the Gaussian 
distribution of the random fields, to show the roughening of the domain interface. Overhangs 
and islands are created. In particular, due to the existence of the intermediate states of the 
spin Si, intermediate regimes emerge. The microscopic structure of the interface becomes 



more complicated than that of the DRFIM model lOl. 111]. 
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There may be different ways to define the domain interface. We adopt the definition 
with the magnetization. Other reasonable definitions yield similar results. Denoting a spin 



at site (x, y 
component 



by S x Jt ) = (cos9 xy (t),smO xy (t)), we introduce a line magnetization of the x 
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m(y,t) 



1 



.x=l 



The height function of the domain interface is then defined as 



Following Refs. 0, H, |ll|, |40|, |41], to describe the propagating and roughening of the domain 
interface, we define the average velocity v (t), roughness function 0J 2 (t), and height correlation 
function C(r,t) of the domain interface respectively, 



h(y,t) = ^-[m(y,t) + l\. 



(2) 



(3) 



v(t) = 



(4) 

(5) 
(6) 



co 2 (t) = (h(y,t) 2 )-(h(y,t)) 2 , 

C(r,t) = ([h(y + r,t)-h(y,t)] 2 ). 

Here (• ■ ■) represents not only the statistic average but also the average in the y direction. 
Further, we introduce the Fourier transform of the height correlation function, 



h(k,t) 

and define the structure factor 



h(y,t)exp(iky), 



(7) 



y y=i 



S(k,t) = (h(k,t)h(-k,t)). 

In Eq. ([8]), (• • •) represents only the statistical average. 

To independently estimate the dynamic exponent z, we introduce a quantity 

F(t) = [M^(t)-M(t) 2 }/io 2 (t). 



(8) 



(9) 



Here M{t) is the global magnetization and M^ 2 \t) is its second moment. In fact, F(t) is 
nothing but the ratio of the planar susceptibility and line susceptibility. 
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B. Scaling analysis 



For the uniform distribution with A > A c and Gaussian distribution of the random 
fields, the depinning phase transition is of second-order. Here A c is the critical strength 
of the random fields separating the first-order phase transition and second-order one. For 
the depinning transition, the order parameter is the average velocity v(t) of the domain 
interface, and it should obey the dynamic scaling theory supported by the renormalization- 



group calculations 18|, [l9|, |42j. In the critical regime, there are two spatial length scales 



in the dynamic relaxation, i.e., the nonequilibrium spatial correlation £(t) ~ t x l z and the 



lattice size L y , and scaling arguments lead to the scaling form [18|, [19|, |42j, 



v{t, r, L v ) = b-? /u G(b- z t, b 1/u r, b^Ly). (10) 

Here, r = (H — H c )/H c , H c is the depinning transition field, b is an arbitrary rescaling 
factor, z is the dynamic exponent, /3 and v are the static exponents. This dynamic scaling 
form is assumed to hold already in the macroscopic short-time regime, after a microscopic 
time scale t m i C , which is sufficiently large in the microscopic sense. Taking b ~ t 1 / 2 , the 
dynamic scaling form is rewritten as 

v(t, r, L y ) = r^ z G(l, t^ z r, r x ' z L y ). (11) 

In the short-time regime, i.e., the regime with £(£) L y , the finite-size effect is negligible, 

v (t,r) =t-P /uz G(t 1/uz r). (12) 

At the transition field H c , i.e., r = 0, a power-law behavior is achieved, 

v(t) ~ t~? /uz . (13) 

With Eq. fll2p . one may locate the transition field H c by searching for the best power-law 
behavior of the interface velocity v(t,r) Q]- To determine v, one derives from Eq. (fl2"|) 

d T ]nv(t,T)\ T=Q ~t 1 ''" 1 . (14) 

In general, even at the transition point, i.e., r = 0, the roughness function u 2 (t) and height 
correlation function C(r,t) may not obey a perfect power- law behavior in early times. Due 
to the random updating scheme in numerical simulations, the domain wall also roughens 
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even without quenched disorder (A = 0, or a = 0). This leads to corrections to scaling. 
To capture the dynamic effects of the disorder, therefore, we introduce the pure roughness 
function and pure height correlation function respectively, 

Du 2 (t) = u 2 (t) -co 2 b (t), (15) 

DC{r,t)=C{r,t)-C b (r,t), (16) 
where u%(t) and C&(r, i) are the roughness function and height correlation function without 



disorder. At the transition point r = and for a sufficient 
DC(r,t) should obey the power-law scaling forms 
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y large lattice L y , Du 2 (t) and 



Duo 2 {t) ~ (17) 

and 

{^2(C-Cioc)/2 r Kio C if r -C £(t) 
* SW . (18) 
t 2 ^ z if £(t) < r 

Here £(t) ~ C is the global roughness exponent, and (i oc is the local roughness exponent. 



The structure factor generally scales as [40j, |4l| , 



S(k,t) = k-& +1 ')f a (kt 1 ' z ), (19) 
where the scaling function takes the form, 

f u 2 (C-<«) if u » 1 

^ , (20) 

[ u 2C+1 if u < 1 

and ( s is the spectral roughness exponent. 

Since oj 2 {t) represents the fluctuation in the x direction and M^(t) — M(t) 2 includes 
those in both the x and y directions, the dynamic scaling behavior of F(t) should be 

F(t) ~ ffl/Ly ~ t^/Ly. (21) 

III. MONTE CARLO SIMULATION 

Our main results are presented for the 6-state DRFCM model, and those for other p are 
also discussed. 
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A. Uniform distribution of random fields 



We first study the depinning phase transition with the uniform distribution of the random 
fields. In Fig. [2(a), the interface velocity v(t) is displayed for different strengths of the 
random fields and different driving fields. For example, at A = 2.5, it drops rapidly for a 
smaller H, while approaches a constant for a larger H. By searching for the best power-law 
behavior, one locates the transition field H c = 0.935(1) and 0.951(1) for A = 1.75 and 2.5, 
respectively. Especially, the microscopic time scale t mic here is very small, about 10 Monte 
Carlo time steps. According to Eq. f fT3|) . one measures the exponent fifvz — 0.251(2) and 
0.208(2) from the slopes of the curves for A = 1.75 and 2.5, respectively. The power-law 
behavior in Fig. [2(a) indicates that the depinning transition is of second order. 

For a smaller A, the velocity approaches a constant for H > A, while drops rapidly to 
zero once H is below A. The power-law behavior in Eq. (TIB"]) could not be detected. This 
is a signal that the depinning transition is of first order, and the transition field H c = A. 
The critical strength A c of the random fields, separating the first-order depinning transition 
and second-order one, varies with the state number p of the spin Si, as shown in Table |T] 
and Fig. [2(b). In fact, this transition at A c could be also detected from the change of the 
morphology of the domain interface 37J, |38|. For the 2D DRFIM model, as A increases, the 



morphology of the domain interface changes from faceted to self-similar at A c = 1. 

The transition field H c varies with p and A. In Fig. [3(a), H c is displayed for different p 
at A = 1.5, and the results are also listed in Table [H For p = 2, i.e., the DRFIM model, 
the transition field H c = 1.2933 agrees with that in Ref. 10j. For a large p, e.g., p = 36, 
H c becomes very small, H c < 0.001. One expects that H c — > in the limit p — > oo. 
This is consistent with the fact that the depinning transition does not exist in the 2D XY 
model with quenched random fields. The p-dependence of H c is different for even and odd 
p because of the asymmetric initial state for odd p. A similar behavior is also observed for 
the p-dependence of A c in Fig. [2(b). We will come back to this point in Subsec. Ill C. In 
Fig. [3(b), the transition fields H c is plotted for different strengths of the random fields, and 
parts of the results are also given in Table [TTJ Different from the case of p = 2, H c of the 
6-state DRFCM model exhibits a non-monotonous and unusual behavior with respect to the 
strength of the random fields. At A = 2.2, H c reaches the maximum. 

In Fig. H(a), the quantity F(t) defined in Eq. ([9]) is shown for A = 1.75 and 2.5 at H c . A 
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power-law behavior is detected but with certain corrections to scaling up to a microscopic 
time scale t m i C ~ 100 or 200, probably due to the multiple states of the spin Si. According 
to Eq. (121 p . a direct measurement from the slopes of the curves give 1/z = 0.743(5) and 
0.764(6) for A = 1.75 and 2.5 respectively. The correction to scaling at early times is 
not universal, but usually in a power-law form. After introducing an empirical power-law 
correction to scaling, i.e., F(t) ~ t l l z {l + c/t), the fitting to the numerical data extends to 
a larger interval, yield 1/z = 0.758(6) and 0.774(6) respectively, and improve the results by 
about two percent. 

To calculate the logarithmic derivative d T lnv(t, r) = d T v(t,r)/v(t,r), we interpolate 
v(t, t) between different H, e.g., in the interval [0.927, 0.943] for A = 1.75, and [0.946, 0.956] 
for A = 2.5. In Fig. 0Jb), d T lnv(t, r) is plotted for A = 1.75 and 2.5 at H c . A power-law 
behavior is observed for t > 100. According to Eq. ( jT4"|) . from the slopes of the curves one 
measures 1/uz = 0.699(9) and 0.683(9) for A = 1.75 and 2.5 respectively. With a power-law 
correction to scaling, the fitting to the numerical data extends to early times, although the 
fitted exponents remain almost unchanged, 1/vz = 0.703(9) and 0.689(9) respectively. 

In Fig. [5]^a), the pure roughness function Du 2 (t) is displayed for different strengths of the 
random fields at H c . A power- law behavior is observed for t > 100. According to Eq. ffTTj) . 
a direct measurement from the slopes of the curves gives 2(/z = 1.74(2) and 1.64(2) for 
A = 1.75 and 2.5 respectively. With a power-law correction to scaling, the results are 
refined to be 2(/z = 1.841(7) and 1.693(7) respectively. 

In Fig. [5]^b), the structure factor S(k, t) as a function of k is displayed for A = 2.5 at H c . 
Obviously, for u = kt 1 ^ <C 1, u^^ +1 ^ f s (u) approaches a constant. For u ^> 1, u^^^ f s (u) 
exhibits a power-law behavior. From the slope of the curve, one extracts the exponent 
2( s + 1 = 2.71(2), i.e., the spectral roughness exponent ( s = 0.86(1), by Eqs. fl]I]) and fl20|) . 
According to Eq. ( fl9|) . data collapse of S(k,t) = S(k,t)t~^ +1 ^ z for different t is shown in 
the inset. The exponents used for the data collapse are £ = 1.09 and z = 1.29 obtained from 
Fig. HUa) and[5]^a) respectively. Since ( s is smaller than the high-/c mode of S(k,t) grows 
with time by a power law, as Eqs. ( fl9l) and (1201) indicate. In other words, the anomalous 
dimension of the fluctuations in a small scale is different from that in a large scale. 

In Fig. [6]^a), the pure height correlation function DC(r,t) is displayed for A = 2.5 at H c . 
For r ^> £(t), e.g., r = 512, one extracts the exponent 2(/z = 1.689(9) by Eq. ffl8|) . consistent 
with that from Fig. [5]^a). For r £(t), e.g., r = 2, DC(r,t) should be independent of t for 
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a normal interface with ( = (i oc , according to Eq. ( TT8|) . But DC(r,t) of r = 2 in Fig. |6]^a) 
clearly increases with time t. A power-law behavior is observed. According to Eq. f fT8|) . the 
slope of the curve gives 2(( — (i oc )/z = 0.699(5). From the measurements of ( and z from 
Fig. H](a) and 0a), we calculate the local roughness exponent (i oc = 0.64(1). Since Q oc is 
smaller than (, DC(r, t) for a small r grows with time by a power law, according to Eq. (1T8|) . 
This behavior is in accord with that of the high-fc mode of S(k,t). We empirically observe 
C ~~ Ooc ~ 2(C — Cs)j although we are not able to theoretically derive it. 

In Fig. Mjo), the pure height correlation function DC(r,t) is plotted as a function of r. 
For r -C £(t), there exists a power-law behavior, but this regime is rather limited. From the 
slope of the curves one estimates 2([ oc = 1.27(2), i.e., (i oc = 0.64(1). For a large r, there 
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emerge strong corrections to scaling. With a specific form of corrections to scaling 

DC(r) ~ [tanh(r/c)] 2Cioc , (22) 

the fitting to the numerical extends to a larger interval, and yields (i oc = 0.639(5), in 
agreement with that from the slope of the curves. This result is also consistent with that 
from Fig. [6]^a) within errors. 

From the measurements of f3/vz, l/z, 1/vz, 2(/z, 2( s + l, 2(£ — Q oc )/z and 2Q oc in Fig. [2] 
to M respectively, we calculate the static critical exponents (3 and v, the dynamic critical 
exponent z, the roughness exponents (, Q oc and ( s for the 6-state DRFCM model. All the 
results are summarized in Table [Til compared with those of the QEW equation. Obviously, 
the critical exponents vary with the strength of the random fields. The critical exponent v 
increases, and all other critical exponents decrease with the strength of the random fields, 
staring from A = 1.25. In general, the critical exponents of the DRFCM model differ from 
those of the QEW equation by 10 to 30 percent. In our opinion, this difference is induced 
by overhangs and islands, and also intermediate regimes created by the intermediate states 
of the spin 5", 1(], 11]. 



At this stage, the dependence of the critical exponents on the strength of the random 
fields for the 6-state DRFCM model is not much different from that for the case of p = 2, i.e., 
the DRFIM model |46l] , although the microscopic structure of the domain interface including 
overhangs and islands, and intermediate regimes, becomes more complicated. This might be 
related to the definition of the height function in Eq ([3]) which averages over the orientations 
of the spins. 
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For the roughening interfaces, according to Ref. |40|, the roughness exponents (, Q oc and ( s 
are not always independent, and there are four different subclasses, namely, Family- Vicsek, 
super-rough, intrinsic and faceted. In Ref. 4jJ, however, it suggests that there is a new 
subclass of anomalous roughening, with ( ^ ( [oc ^ ( s and (i oc ^ 1. Our results in Table HT1 
indicate that the growth dynamics of the domain interface in the DRFCM model, similar to 
that in the DRFIM model [46|, belongs to such a subclass. 

Looking at Fig. [TJ one might catch such a scenario that the domain- wall motion could be 
decomposed into modes: one is the growth of the width of the domain wall, and another is 
the roughening of the domain wall. The subtle point here is how to define the width of the 
domain wall. Essentially, the width of the domain wall should be generated by overhangs 
and islands, and intermediate regimes. Therefore, one possibility is to define the width of the 
domain wall by the size of overhangs and islands, and intermediate regimes. For p = 2, there 
are no intermediate regimes, and the dynamic relaxation of overhangs and islands have been 
carefully analyzed in Ref. The average size of overhangs and islands does grow with 

time by a power law. Denoting the velocity of overhangs and islands by 5v , and the velocity 
of the domain wall defined by the magnetization by Vm, one obtains Sv/vm ~ t e , with 
6 ~ 0.5 [11]. For p=6, one may similarly define the width of the domain wall, although due 
to the existence of intermediate regimes, the fluctuation is somewhat larger. The exponent 
9 is estimated to be 6 « 0.25. Since the paper is already lengthy, details of these contents 
will not be presented here. 

However, it is still not possible to find a decomposition scheme such that the roughening 
mode of the domain wall is the same as that of the QEW equation. In some sense, the 
height function defined by the magnetization already minimizes the effect of overhangs and 
islands, and intermediate state, but its dynamic evolution is still affected by overhangs and 
islands, an intermediate states. 



B. Gaussian distribution of random fields 



We now consider the Gaussian distribution of the random fields. As shown in Fig. [3] and 
Tables HI and IHIl the transition field H c varies also with the state number p of the spin Si and 
the standard deviation a of the Gaussian distribution. Since the Gaussian distribution is 
unbounded, the crossover of the depinning transition from the second order to the first order 
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is not present. Qualitatively, the functions H c (p) and H c (a) for the Gaussian distribution are 
similar to those for the uniform distribution. In particular, if one makes a correspondence 
A — > 1.5a, the transition field H c for the Gaussian distribution looks rather close to that 
for the uniform distribution, as shown in Fig. [3] and Tables HI and IIHI Quantitatively, H c (a) 
for the Gaussian distribution fits nicely to a parabolic curve as can be seen in Fig. E^b). 

The critical exponents for the 6-state DRFCM model with the Gaussian distribution of the 
random fields are given in Table [TTT1 Although the dependence of the critical exponents on the 
strength of the random fields looks similar for both the Gaussian and uniform distributions 
in the regimes A > 1.5 and o > 1, the critical exponents do vary with the form and strength 
of the random fields. The strong universality is violated. For p = 2, i.e., the DRFIM model, 
some critical exponents or their combinations, e.g., /3 and (2( + l)/z, are 'quasi universal'. In 
other words, these critical exponents are independent of the form and strength of the random 
fields in certain regimes, e.g., A > 1.5 and a > 1. However, such a quasi universality does 
not exist for the 6-state DRFCM model. This should be due to the multiple orientations of 
the spin Si. The dynamic evolution of the domain interface is more sensitive to the quenched 
disorder. 

C. Orientations of ordered domains 

For p = 2, i.e., the DRFIM model, there are only two possible orientations of the mag- 
netization of the two ordered domains, and the domain interface is unique. For p > 4, 
different orientations of the two ordered domains may lead to different domain interfaces. 
For p = 3, there exist two domain interfaces, but they are equivalent. Now, an interesting 
question is, whether the depinning phase transitions of the domain interfaces with different 
orientations of the two ordered domains occur at a same transition field, and belong to a 
same universality class, or not. 

To determine the transition field H c , we again make use of Eq. f jT2|) . In Fig. [7(a), one 
locates H c for the 4-state DRFCM model with the Gaussian distribution of the random fields 
at a — 1.5, by searching for the best power-law behavior of the interface velocity v(t,r). 
Obviously, the transition field H c depends on the orientations of the two ordered domains. 
For 9 L = and 9 R = vr, H c = 1.2833(3); for 9 L = and 6 R = tt/2, H c = 1.5645(3). In 
Table IIV} the results for p = 4, 6 and 8 are summarized. 
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Since the driving field H is applied in the x direction, one may naively think that the 
smaller Or is, the easier it is to drive the domain wall to propagate. Therefore, the transition 
field H c may be expected to decrease with Or. Looking at Table IIV} however, one finds 
rather unusual dependence of H c on Or. As Or decreases, except for Or = n, H c does not 
decreases, rather increases! For p = 8, for example, H c = 0.6018(5) for Or = 37r/4, while 
H c = 1.7825(5) for Or = 7r/4. The difference is very prominent. It is much more difficult to 
drive the domain wall of Or = tt/A to propagate. 

Why does such an unusual phenomenon occur? For Or ^ it, due to the action of the 
driving field H and the average field of the nearest neighbor spins Sj, a spin Si on the 
domain interface rotates almost exclusively to the direction 0^ < Or. In Fig. H](b), for 
example, the snapshot for Or = 2ir/3 is displayed for the 6-state DRFCM model. Almost 
no spins orient in the directions 0i — 7r and 47r/3, and only a few spins in the direction 
Qi = 57r/3. In other words, for Or = 27r/3, spins may be flipped to 0i = tt/3 and 0, while 
for Or = 7r/3, to 0i = only. For Or = tt/3, 0{ is pinned if the flip to 0i = is forbidden 
by the random fields. For Or = 2tt/3, however, there is still another channel to Qi = 7r/3. 
Therefore, H c of the domain interface with Or = tt/3 is larger. 

The transition field H c for Or = tt is smaller than that for the smallest Or, while larger 
than others. One may understand the former phenomenon following the argument above, 
but the latter is puzzling. A possible reason might be that in the case of Or = tt, the average 
field of the nearest neighbor spins Sj in the domain interface is zero, only the external field 
drives the spins to the x direction, therefore, H c is somewhat larger. 

Another issue one still needs to address is that according to Table IIVI the transition 
field H c for the smallest Or increases with p. This seems contradict with the fact that the 
depinning transition is absent in the limit p — > oo. Therefore we compute H c for the smallest 
Or up to p = 24, and the results show that H c starts dropping around p = 16. 

With the procedure in Subsec. Ill A, we may determine all the critical exponents for 
the domain interfaces with different orientations of the two ordered domains. In Fig. [7(a) 
and (b), for example, one measures the exponent (3/vz from the slopes of the curves v(t) 
at H c . From the results summarized in Table IIV| we observe that the critical exponents 
do change with the orientations of the two ordered domains, and this phenomenon tends 
to be more prominent for larger p. In other words, the depinning phase transitions of the 
domain interfaces with different orientations of the two ordered domains belong to different 
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universality classes. As Or decreases, the dynamic exponent z, roughness exponents (, 
(ioc and ( s decrease monotonously, while the static exponent v increases. The dependence 
of the static exponent /3 on Or is somewhat complicated. Roughly speaking, a larger H c 
corresponds to a smaller /3. The exponent /3/v is, however, almost independent of Or, except 
for the smallest Or. 

Comparing the results in Table IIVI with those in Tables [TT] and 1111} one may find that as 
Or decreases, the dynamic effect of the random fields tends to be stronger. In other words, 
for smaller Or, 0{ is more sensitive to the random fields. In any case, the domain interfaces 
always belong to the new subclass of anomalous roughening, with ( ^ Q oc ^ ( s and Q oc ^ 1. 



IV. CONCLUSION 



With Monte Carlo simulations, we have investigated the depinning phase transition of 
the domain interface in the 2D p-state DRFCM model, emphasizing the dynamic effect 
of the orientations of the magnetization of the two ordered domains. Based on the short- 
time dynamic approach, the transition field and all the critical exponents are determined. 
Numerical results are mainly presented for p = 6, and also extended from p = 2 to 12. 

The important results are summarized below: 

i) For the uniform distribution of the random fields, there exists a critical strength A c of 
the random fields, separating the first-order depinning transition (A < A c ) and second-order 
one (A > A c ). For the Gaussian distribution, the depinning transition is of second order. 
For a fixed strength A of the random fields, the transition field H c — > in the limit p — > oo. 

ii) For the second-order depinning transition, the critical exponents do vary with the 
form of the random-field distribution and the strength of the random fields, and the strong 
universality is violated. Further, we observe ( ^ Q oc ^ ( s , and Q oc < 1, and it indicates that 
the r oug hening dynamics of the domain interface belongs to the new subclass suggested in 



Ref. 



4l|. 



iii) More importantly, the transition field H c and critical exponents change with the initial 
orientations of the magnetization of the two ordered domains. H c exhibits rather unusual 
dependence on the orientation Or as shown in Table IIVI Unexpectedly, the dynamic effect 
of the random fields is the most significant for the domain interface with the smallest Or. 
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p 


A c 


H c 


P 


A c 








A = 1.5 cr = 1.0 






A = 1.5 a = 1.0 


2 


1.000(1) 


1.2933(2) 1.4205(3) 


8 


0.42(1) 


0.757(1) 0.778(2) 


3 


1.000(1) 


1.341(1) 1.513(2) 


9 


0.53(1) 


0.594(2) 0.607(2) 


4 


1.000(1) 


1.192(1) 1.255(1) 


10 


0.41(1) 


0.634(2) 0.664(2) 


5 


0.625(5) 


0.906(1) 0.968(2) 


11 




0.515(2) 0.526(2) 


6 


0.51(1) 


0.921(1) 0.946(1) 


12 




0.553(3) 0.585(3) 


7 


0.58(1) 


0.705(1) 0.727(2) 


36 




< 0.001 < 0.001 



TABLE I: The critical strength A c for the uniform distribution of the random fields, and the 
depinning transition field H c at A = 1.5 and a = 1.0 for the uniform and Gaussian distributions 
of the random fields respectively. 





v(t) 


F(t),C(r,t) 


u 2 (t),C(r,t) 


C(r,t) 


S(k,t) 


A 


H c 


P 


;/ 


P/v 


z 


c 


Cloc 


Cs 


1.25 


0.913(1) 


0.491(8) 


0.95(2) 


0.517(4) 


1.40(1) 


1.24(1) 


0.78(1) 


1.05(1) 


1.5 


0.921(1) 


0.415(8) 


1.04(2) 


0.399(4) 


1.37(1) 


1.22(1) 


0.75(1) 


1.00(1) 


1.75 


0.935(1) 


0.358(6) 


1.06(2) 


0.338(3) 


1.33(1) 


1.21(1) 


0.72(1) 


0.97(1) 


2.0 


0.954(1) 


0.333(6) 


1.08(2) 


0.308(3) 


1.32(1) 


1.18(1) 


0.70(1) 


0.92(1) 


2.25 


0.963(1) 


0.313(6) 


1.10(2) 


0.285(3) 


1.31(1) 


1.13(1) 


0.67(1) 


0.89(1) 


2.5 


0.951(1) 


0.302(6) 


1.12(2) 


0.270(3) 


1.29(1) 


1.09(1) 


0.64(1) 


0.86(1) 


DRFIM 


1.2933(2) 0.295(6) 


1.02(2) 


0.289(3) 


1.33(1) 


1.14(1) 


0.74(1) 


0.94(1) 


QEW 




0.33(2) 


1.33(4) 


0.25(2) 


1.50(3) 


1.25(1) 


0.98(6) 


1.25(1) 



TABLE II: The depinning transition field and critical exponents of the 6-state DRFCM model for 
different strengths of the random fields with the uniform distribution are compared with those of 



the QEW equation [4 



and DRFIM model with A 



1.5 QQ. 
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V (t) 


F(t),C(r,t) 


co 2 (t),C(r,t) 


C(r,t) 


S{k,t) 


a 


H c p » P/v 


z 


C 


Cloc 


Cs 


1.0 


0.946(1) 0.380(6) 1.00(2) 0.380(3) 


1.35(1) 


1.26(1) 


0.76(1) 


0.97(1) 


1.2 


0.965(1) 0.337(6) 1.05(2) 0.321(3) 


1.33(1) 


1.19(1) 


0.73(1) 


0.90(1) 


1.4 


0.968(1) 0.321(6) 1.13(2) 0.284(3) 


1.31(1) 


1.12(1) 


0.69(1) 


0.83(1) 


1.6 


0.957(1) 0.312(6) 1.18(2) 0.264(3) 


1.29(1) 


1.06(1) 


0.65(1) 


0.80(1) 



TABLE III: The depinning transition field and critical exponents of the 6-state DRFCM model for 
different strengths of the random fields with the Gaussian distribution. 







v(t) 


F(t),C(r,t) 




CM) 


S(k,t) 




L Or 


H c 


P 


V 


P/v 


z 


C 


Cloc 


Cs 


p = 4 





7T 


1.2833(3) 


0.319(6) 


1.10(2) 


0.290(3) 


1.32(1) 


1.11(1) 


0.63(1) 


0.81(1) 









1.5645(3) 


0.280(6) 


1.16(2) 


0.241(3) 


1.26(1) 


1.05(1) 


0.60(1) 


0.76(1) 


p = 6 





7T 


0.9643(5) 


0.319(6) 


1.15(2) 


0.277(3) 


1.30(1) 


1.10(1) 


0.66(1) 


0.81(1) 







h 


0.8020(5) 


0.325(6) 


1.17(2) 


0.278(3) 


1.26(1) 


1.03(1) 


0.57(1) 


0.72(1) 







h 


1.6808(5) 


0.287(6) 


1.18(2) 


0.243(3) 


1.25(1) 


1.02(1) 


0.56(1) 


0.67(1) 


p = 8 





TT 


0.7841(5) 


0.314(6) 


1.15(2) 


0.273(3) 


1.39(1) 


1.15(1) 


0.65(1) 


0.90(1) 







h 


0.6018(5) 


0.322(6) 


1.16(2) 


0.278(3) 


1.34(1) 


1.08(1) 


0.57(1) 


0.83(1) 







¥ 


0.6895(5) 


0.334(6) 


1.21(2) 


0.276(3) 


1.26(1) 


1.01(1) 


0.54(1) 


0.65(1) 









1.7825(5) 


0.308(6) 


1.25(2) 


0.246(3) 


1.21(1) 


0.98(1) 


0.52(1) 


0.63(1) 



TABLE IV: The depinning transition field and critical exponents of the p-state clock model with 
the Gaussian distribution at a = 1.5 for different initial orientations of the magnetization of the 
two ordered domains. 
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(a) cr = 1.5, if = 0.9643 



(b) cr = 1.5, H = 0.8020 



FIG. 1: (Color online) Snapshots of the spin configurations at t = 500 Monte Carlo time steps, for 
the 6-state DRFCM model with the Gaussian distribution of the random fields at H r . Arrowheads 



denote the orientation of the spins, (a) 6l = and Or = 7r; (b) 9l = and Or = 2tt/3. 




FIG. 2: (Color online)(a) The interface velocity v(t) is plotted for A = 2.5 at different driving field 
H with solid lines. For comparison, v(t) is also shown for A = 1.75 and 2.5 at H c = 0.935 and 
0.951 respectively with circles and crosses. Dashed lines show power-law fits. For clarity, the solid 
lines are shifted up by a factor of five, (b) The critical strength A c of the disorder is displayed for 
different p. In both (a) and (b), the results are for the uniform distribution of the random fields. 
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FIG. 3: (Color online) (a) The depinning transition field H c is displayed for different p at A = 1.5 
and a = 1.0 for the uniform and Gaussian distributions of the random fields respectively. Errors 
are smaller than the symbol sizes, (b) H c is plotted for different strengths of the random fields. In 
the inset, H c for p = 2, i.e., the DRFIM model, is shown for comparison. 




FIG. 4: (Color online) (a) The quantity F(t) is shown for A = 1.75 and 2.5 at H c . (b) The 
logarithmic derivative of v(t,r) is plotted for A = 1.75 and 2.5 at H c . In both (a) and (b), the 
results are for the uniform distribution of the random fields. Dashed lines show power-law fits, and 
solid lines include corrections to scaling. For clarity, the curves for A = 2.5 are shifted up by a 
factor of three. 
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FIG. 5: (Color online) (a) The pure roughness function Dui 2 (t) is displayed for A = 1.75 and 2.5 
at H c . Solid lines represent power-law fits with corrections to scaling. For clarity, the curve for 
A = 2.5 is shifted up by a factor of five, (b) The structure factor S(k,t) is plotted for A = 2.5 at 
H c . In the inset, data collapse for S(k,t) = S(k,t)t~( 2< * +1 ^ z for different t is shown. In both (a) 
and (b), the results are for the uniform distribution of the random fields, and dashed lines show 
power-law fits. 
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FIG. 6: (Color online) (a) The pure height correlation function DC(r,t) is displayed for A = 2.5 



at H c . The solid line shows a powerdaw fit with a correction to scaling, (b) DC(r,t) is plotted as 
a function of r for A = 2.5 at H c . The solid line shows a powerdaw fit with a hyperbolic tangent 
correction. In both (a) and (b), the results are for the uniform distribution of the random fields, 
and dashed lines show powerdaw fits. 
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FIG. 7: (Color online) (a) The interface velocity v(t) is plotted for the 4-state DRFCM model with 
the Gaussian distribution of the random fields at different driving field H with solid lines, for the 
initial orientation 9ft = tt. For comparison, v(t) is also shown for Or = tt and tt/2 at H c respectively 
with crosses and circles. For clarity, the solid lines are shifted up by a factor of three, (b) v(t) is 
displayed for the 6-state DRFCM model at H c , with different Or. In both (a) and (b), dashed lines 
represent power-law fits. 
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